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We report on several improvements of the local bosonic algorithm proposed by M. Liischer. We find that 
preconditioning and over-relaxation works very well. A detailed comparison between the bosonic and the Kramers- 
algorithms shows comparable performance for the physical situation examined. 
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1. INTRODUCTION 

Since its proposal in M, the local bosonic 
fcrmion method has been thoroughly investigated 
and improved §-|,§-|l|. We report here on 
some improvements, especially the efficient im- 
plementation of preconditioning techniques, and 
give more detailed information on the comparison 
of this algorithm to the Kramers algorithm pre- 
sented in [fl5| . As in ||, all investigations were 
made with the gauge group SU(2). 

2. PRECONDITIONING 

Even-odd preconditioning is known to work 
very well for iterative solvers. It is therefore 
straightforward to try to apply this technique to 
the local bosonic algorithm. We quickly recall 
that the hermitean dirac operator can be written 
as 
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Co = [ca/(1 + 8K)]~ 1 with cm > 1 is a normaliza- 
tion constant so that \\Q\\ < 1. Using the follow- 
ing identity, 
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we can immediately write down the precondi- 
tioned matrix 
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with Cq = [cm(1 + MK 2 )}- 1 and ||Q|| < 1. This 
matrix, however, contains next-to-nearest neigh- 
bour interactions. Since it appears squared in the 



local bosonic action, the local updates become 
rather complicated. The local bosonic approxi- 
mation to the fermion determinant is 
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with P a suitably chosen approximation to 1/x. 
By applying Equation (^) to each single factor in 
this equation, we get 
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Letting Q = Co/cqQ, we obtain the precondi- 
tioned action 



S b = J24(Q ~ Porl)(Q ~ P r k ) 



(6) 



which is very similar to the original one. P de- 
notes the projector on the odd sites. 

The dynamics of the preconditioned algorithm 
is investigated in Table [l| As can be seen, pre- 
conditioning saves a factor of about 4 in computer 
time, as expected. The analysis of the condition 
numbers of Q 2 and Q 2 show that the gain is even 
near 8. The reason for that is that the largest 
eigenvalue of Q 2 is far below 1. The attempt 
to exploit this fact by adjusting cm fails because 
when the largest eigenvalue of Q 2 gets near 1 the 
bosonic fields develop slow modes. Detailed in- 
vestigations show that there is a tradeoff between 
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Table 1 

Autocorrelation times in units of 1OOOQ0- multiplications on the 4 3 x 8 lattice for (3 = 1.75 and K — 0.165. 
n and e are chosen so that 6 = 0.03. All except the first run are preconditioned. The letters in the third 
column give the type and order of sweeps used per iteration, where H is a bosonic heatbath, O a bosonic 
over-relaxation and h and o the gauge updates. 



the slow bosonic modes and the number of fields. 
This is however not yet analytically understood. 

The idea of applying some reverse transforma- 
tion to the single factors of the polynomial can 
easily be applied to different cases of precondi- 
tioning, e.g. preconditioning with the inverse 
free fermion matrix, still resulting in a local al- 
gorithm. 

3. HYBRID OVER-RELAXATION 

It is believed that the application of local heat- 
bath and over-relaxation sweeps in a combination 
1 : £, where £ is a suitably defined correlation 
length, results in a dynamical critical exponent 
of z w 1 (see e.g. 0). It was found in || that 
the application of this method to the gauge fields 
or the bosonic fields only has no effect, which can 
be understood as the fact that the coupling of the 
fields dominates the dynamics of the algorithm in 
the examined case. The application to both fields 
however is very efficient, as can be seen from Ta- 
ble ||. Although this is far from a verification of 
2wl, we believe that the gain will increase with 
increasing correlation length. 

4. COMBINED UPDATE 

In |s| an update which updates the bosonic and 
gauge fields at the same time was proposed. The 
method uses the fact that we can update a gauge 
link U(x,fi) according to the effective action af- 
ter integrating out all bosonic degrees of freedom 
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Table 2 

Autocorrelation times in units of lQ0QQcf>- 
multiplications. The run parameters are the same 
as in the third line of Table 1. 



4>(x) and (j){x + (x) and then refresh </>(x) and 
<p(x+p,) with heatbath-updates. The hope is, that 
the effective action induces a smaller fermionic 
force on the gauge field, which was found to dom- 
inate the update. Indeed, it is only of order K 2 
instead of K in the standard update; however 
the updating procedure is more complicated what 
may compensate for this effect. There also exists 
no free-field analysis to get a hint about the dy- 
namical critical exponent of this update. Table 
H shows that the update performs comparably to 
the standard updates. However its implementa- 
tion is much more difficult and its application to 
improved actions is probably not feasible. 

5. ALGORITHM COMPARISON 

In JOI the comparison between the local 
bosonic algorithm in its hcrmitcan version and 
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Table 3 

Autocorrelation times for the combined update. 
The first line uses one combined over-relaxation 
step, for the second a standard bosonic over- 
relaxation step was added. 
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Table 5 
Results for both algorithms. 



the Kramers algorithm has been described. We 
will give more detailed information about the 
runs made and the parameters chosen. The 
bosonic parameters were chosen according to the 
investigations presented above while the Kramers 
algorithm was examined in Q . In Table ^ the pa- 
rameters for both algorithms are given as well as 
the expectation values of the extremal eigenvalues 
of Q 2 (with cm given as in the Boson section of 
the table) . In Table || the results for the observ- 
ables are presented. As can be seen, the param- 
eters were chosen in a way that the observables 
agree; it seems however that 8 is slightly too big 
in the 8 3 12 case. 
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